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Abstract 

In the framework of braneworld formalism, we examine the dynamics of a 
Bianchi IX model with three scale factors on a 4-dim Lorentzian brane embed¬ 
ded in a 5-dim conformally flat empty bulk with a timelike extra dimension. The 
matter content is a pressureless perfect fluid restricted to the brane. In this sce¬ 
nario Einstein’s equations on the brane reduces to a 6-dim Hamiltonian dynamical 
system with additional terms - due to the bulk-brane interaction - that avoid 
the singularity and implement nonsingular bounces in the early phase of the uni¬ 
verse. Due to an effective cosmological constant on the brane the phase space of 
the model presents two critical points (a saddle-center-center and a center-center- 
center) in a finite region of phase space, and two asymptotic de Sitter critical 
points at infinity, one acting as an attractor to late-time acceleration dynamics. 
The critical points belong to a 2-dim invariant plane; together they organize the 
dynamics of the phase space. The center-center-center critical point corresponds 
to a stable Einstein universe configuration with perpetually oscillatory orbits in 
its neighborhood. On the other hand the saddle-center-center engenders in the 
phase space the topology of stable and unstable 4-dim cylinders R x S' 3 , where 
R is a saddle direction and S' 3 is the center manifold of unstable periodic orbits, 
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the latter being the nonlinear extension of the center-center sector. By a proper 
canonical transformation we are able to separate the three degrees of freedom of 
the dynamics into one degree connected with the expansion and/or contraction of 
the scales of the model, and two pure rotational degrees of freedom associated with 
the center manifold S 3 . It follows that the typical dynamical flow is an oscillatory 
mode about the orbits of the invariant plane. For the stable and unstable cylinders 
we have the oscillatory motion about the separatrix towards the bounce, leading 
to the homoclinic transversal intersection of the cylinders, as shown numerically in 
two distinct experiments. We show that the homoclinic intersection manifold has 
the topology of R x S 2 consisting of homoclinic orbits biasymptotic to the center 
manifold S 3 . This behavior defines a chaotic saddle associated with S 3 , indicating 
that the intersection points of the cylinders have the nature of a Cantor set with 
a compact support S' 2 . This is an invariant signature of chaos in the model. We 
discuss the possible connection between these properties of the dynamics, namely 
the oscillatory approach to the bounce together with its chaotic behavior, and 
analogous features present in the BKL conjecture in general relativity. 

PACS numbers: 05.45.-a, 05.45.Pq, 98.80.-k, ll.25.-w 


1 Introduction 

The general Bianchi IX model has become a paradigm for the behavior of 
the general relativity dynamics near the cosmological singularity since the 
seminal papers of Belinskii, Khalatnikov and Lifshitz (BKL)[1 ) j2[ 3j, and 
collaborators[U [5]. They showed that in a Bianchi IX model with three scale 
factors the approach to the singularity (t —> 0) is an oscillatory mode, con¬ 
sisting of an infinite sequence of Kasner eras in each of which two of the 
scale factors oscillate while the third decreases monotonically. On passing 
from one era to another (with decreasing time t ) the monotonic behavior 
is transferred to another of the three scale factors. It was also shown that 
(i) the length of each era was determined by a sequence of numbers x s , 
0 < x s < 1, s = integer, each of which arises from the preceding one by the 
map z s _(_i = the fractional part of l/x s , with the length of the s-th era given 
by k s = the integral part of l/x s ; (ii) this map leads to spontaneous stochas- 
tization in the sequence of eras on approaching the singularity (t —> 0) for 
arbitrary initial conditions given at t > 0. Due to the involved nonintegrable 
dynamics the evolution of the model had to be actually treated in asymptotic 
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regions of arbitrarily small times together with truncations made to guaran¬ 
tee the validity of the perturbation method, so that “in the most general 
case all details of such regime are not yet fully understood” [Tj . In the past 
four decades the dynamics of these models has been reexamined in an exten¬ 
sive literature but - as in the BKL work - the approach has been basically 
twofold: to obtain maps which approximate the dynamics and which exhibit 
strong stochastic properties, and the discussion of how well these discrete 
maps represent the full nonlinear dynamics of Bianchi IX models in general 
relativity. From the point of view of the phase space flow the interest in the 
chaoticity of Bianchi IX models has been mainly focused on the mixmaster 
universe (the vacuum Bianchi IX case with three scale factors), although the 
question of the behavior (chaotic or not) remained unsettled mainly due to 
the absence of an invariant characterization of chaos in the model (standard 
chaotic indicators as Liapunov exponents being coordinate dependent and 
hence questionable). Therefore, along with the Cosmic Censorship Conjec¬ 
ture, the BKL conjecture is probably one of the major unsolved issues of 
classical general relativity connected to the presence of a singularity in the 
dynamics. 

Our purpose in the present paper is to examine the dynamics of a 4-dim 
Bianchi IX model with three scale factors in the framework of a braneworld 
formalism (which encompasses general relativity as a classical low-energy 
limit). Due to an extra timelike dimension brane-bulk interaction terms 
correct general relativity substituting the singularity by nonsingular bounces 
in the cosmological dynamics. The dynamics of the approach to the bounces 
is extremely complex presenting oscillatory and chaotic features of the BKL- 
type but, as we will show, they are amenable to an exact analytical/numerical 
treatment so that we may have a more clear picture of what happens in the 
general relativity limit. 

Most of the approaches to the problem of the initial singularity and to 
the possible solutions adopted to circumvent this problem lie in the realm 
of a quantum theory of gravitation. In fact we may consider that the initial 
conditions of our present expanding Universe were fixed when the early Uni- 
verse emerged from a Planckian regime and started its classical evolution. 
However, by evolving back the initial conditions using Einstein’s classical 
equations the Universe is driven toward a singular point where the classical 
regime is no longer valid[6|. This is an indication that classical general rela¬ 
tivity is not a complete theory and in this domain quantum processes must 
be taken into account. 
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Among several propositions to describe the dynamics in the semiclas- 
sical domain prior to the classical regime are, for instance, quantum loop 
cosmology [7] and the string based formalism of D-branes[8], both of them 
leading to corrections in Einstein’s equations and encompassing general rel¬ 
ativity as a classical (low energy) limit. In the present paper we adhere to 
the so-called braneworld scenario [5] -[TU] based on the string formalism of D- 
branes. In this context, extra dimensions are introduced by a bulk space and 
all the matter in the Universe would be trapped on a brane embedded in the 
bulk with three spatial dimensions; only gravitons would be allowed to leave 
the brane and move into the full bulk [11] . At low energies general relativity 
is recovered but at high energies significant changes are introduced in the 
gravitational dynamics. Our interest in this framework comes from the fact 
that it can provide corrections that are dominant in the neighborhood of the 
singularity, resulting in a repulsive force which avoids it completely and leads 
the Universe to undergo nonsingular bounces. Bouncing brane world models 
were constructed by Shtanov and Salmi [TO] based upon a Randall-Sundrum 
type action with one extra timelike dimension. A complete analysis of bounc¬ 
ing brane world dynamics embedded in a five-dimensional de Sitter spacetime 
may be found in Refs. pans], where both high energy local corrections as 
well as nonlocal bulk corrections are analyzed on a spatially homogeneous 
brane. 

Although spacelike extra dimensions theories have received more atten¬ 
tion in the last decades [IT], studies regarding extra timelike dimensions 
have been considered [HI [T5l. FT6I T7| . Albeit presenting some problematic 
issues [181 19] it has been shown[20] that they might be circumvented by con¬ 
sidering a noncompact timelike extra dimension, which is the case of the 
model in this paper. 

In our braneworld scenario we consider a 5-dim de Sitter bulk space with 
a timelike extra dimension, and a 4-dim Lorentzian brane with a Bianchi 
IX geometry with three scale factors. The matter content of the model is 
taken as a pressureless perfect fluid (dust) restricted to the brane and an 
effective nonvanishing cosmological constant is also considered. With the 
above assumptions we show that the Gauss-Codazzi equations, and hence, 
the modified field equations on the brane are automatically satisfied. The 
modified Einstein’s equations for the model have a first integral that can be 
expressed as a Hamiltonian constraint, yielding a three degrees of freedom 
dynamical system in a 6-dim phase space. The additional correction terms 
due to the bulk-brane interaction avoid the initial singularity resulting in- 
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stead to nonsingular bounces in the model. One of the main features of the 
phase space is the presence of a saddle-center-center critical point with an 
associated center manifold of unstable periodic orbits having the topology 
S 3 . We will show that from the center manifold it emerges stable and un¬ 
stable manifolds with the topology of spherical cylinders Rx S 3 (constituted 
actually of bounded oscillatory orbits) which cross each other transversally 
in the neighborhood of the bounces. These transversal crossings provide an 
invariant characterization of homoclinic chaos in the model. 

These results are in realm of recent studies in the characterization of 
homoclinic chaos for Hamiltonian dynamical systems with n > 2 degrees of 
freedom. For n = 2 the characterization of chaos connected with the presence 
of homoclinic phenomena in the dynamics has been the object of an exten¬ 
sive outstanding literature (cf. [21] (22, 23jj2Tll25, 26] and references therein). 
The dynamics near homoclinic orbits is very complex, with the homoclinic 
intersection manifold associated with the presence of the well-known horse¬ 
shoe structures (cf. |27[ [28l [29] and references therein), which is an invariant 
signature of chaos. Furthermore, invariant Cantor sets associated with a 
horseshoe construction are connected to chaotic saddles [3D]-[32]. For n > 2 
orbits homoclinic to the center manifold are expected to exist. It has been 
shown, for instance, that critical points of the type saddle-center-...-center 
induce reaction type dynamics in the framework of Transition State Theory 
(see [33] and references therein). The existence of such homoclinic orbits 
has been studied in [33] El] • Although there are no theorems describing the 
dynamics connected with orbits homoclinic to S 3 , it has been shown[35] that 
if there is a transversal intersection of the stable and unstable manifolds, a 
chaotic saddle, and hence an homoclinic trajectory must exist. An inter¬ 
esting analysis of this feature was given in [31], where the authors provide 
a computational procedure to detect a chaotic saddle (and thus homoclinic 
orbits) in the case of Hamiltonian systems with three degrees of freedom. In 
the present paper we follow an alternative procedure to show the presence of 
homoclinic connections with the center manifold S 3 . 

We organize the paper as follows. In the next section we present a brief 
introduction to BraneWorld Theory, deriving the modified field equations on 
the brane. In Section III we construct a general Bianchi IX cosmological 
brane model, with an effective cosmological constant and the matter content 
being dust. In Section IV we study the structure of the phase space, identify¬ 
ing the constants for the linearized motion. In Section V the dynamics about 
the saddle-center-center critical point is examined. Section VI is devoted to 
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a complete analysis of the nonlinear center manifold, together with the 4-dim 
stable and unstable cylinders that emanate from it. Finally in Section VII we 
study the homoclinic transversal intersections of the cylinders that gives an 
invariant characterization of chaos in the dynamics. Conclusions and future 
perspectives are presented in the final section. 

2 The field equations 

For sake of completeness we give here a brief introduction to Braneworld 
Theory, making explicit the specific assumptions used to obtain the dynamics 
of the model. We refer to mm for a more complete and detailed discussion 
and our notation closely follows [6]. We start with a 4-D Lorentzian brane £ 
with metric ^ g ab , embedded in a 5-D conformally flat bulk M with metric 
5 9ab■ Capital Latin indices run from 0 to 4, small Latin indices run from 
0 to 3. We regard £ as a common boundary of two pieces Mi and M 2 of 
Ai and ^g ab is the induced geometry on the brane by the metrics of the two 
pieces. These metrics should coincide on £ although the extrinsic curvatures 
of £ with respect to A4i and M .2 can be different. The action for the theory 
has the general form 

S = 7^2 { f y/-e ®g [ (5) R - 2 A 5 + 2 k 2 5 L 5 \ d 5 x 
Um! 

+ [ V-e ^g [ {5) R - 2A 5 + 2 k 2 5 L 5 ] d 5 x 

Jm 2 

+ 2 ej yj—A)gK 2 d 4 x - 2e J ^/-AgK^xX (1) 
V i Mh (^|‘ 'r - 2 -) ^ 

+ / V-^gL 4 ( {4) g ab: p)d 4 x . 

J s 

In the previous equation, ^ R is the Ricci scalar of the Lorentzian 5-D metric 
^9ab on 2 W, and is the scalar curvature of the induced metric ^g ab on 
£. The parameter cr denotes the brane tension. The unit vector n A is normal 
to the boundary £ and has norm e. If e = +1 the signature of the bulk space 
is so that the extra dimension is timelike. The quantity 
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K = Kg^h/'h is the trace of the symmetric tensor of extrinsic curvature 
K ab = Yf V c n D , where Y A (x a ) are the embedding functions of X in 
While L 4 (^ g ab , p) represents the Lagrangian density of the perfect 
fluid [37] (with equation of state p = ap ), whose dynamics is restricted to the 
brane X, L b denotes the Lagrangian of matter in the bulk. All integrations 
over the bulk and the brane are taken with the natural volume elements 
y—7 g d 5 x and \J—d A x respectively. Einstein constants in five- and 

four-dimensions are indicated with k§ and k\ = 8ttGn, respectively (GV 
being the Newton’s constant on the brane). Throughout this section we use 
natural units with H — c—1. 

Variations that leave the induced metric on X intact, furnish the equations 

® Gab + A 5 ^gAB = ^5 ^ Tab ■ (2) 

Considering arbitrary variations of ® g A s and taking into account Eq. ([2]), 
we obtain 


(i) G ab + e^f ( 5 ^ + S^A = {r ab - ogab) , (3) 

where S ab = K ab — K^g ab , and T ab is the energy momentum tensor on the 
brane. In the limit ^ Kg, Eq. ([3]) reduces to the Israel-Darmois junction 
conditions [38] 


+ S ab ) = eK l i T ab - V (4) 9ab ) ■ ( 4 ) 

Imposing the Z 2 -symmetry[II] and using the junction conditions (Eq. [I]), 
we determine the extrinsic curvature on the brane, 





Now using Gauss equation 



(5) 


(4) iW = (5) R.MNRsY™Y»Y*Y s d - 
e(K m K bd - K aJ K lc ) . 

together with Eqs. ([2]) and ((5]) we obtain the induced field equations on the 
brane 


(4) G ab + A4 {4) 5' afe = 8nG N T ab - eKgII ab + 

tG ab T aF ab . 


(7) 


7 






In the above E ab = ^Cabcd^Y^uFY^ is the projection of the 5-D Weyl 
tensor, and we have defined 


A 4 

Gn 

n a & 


F a b 




HZ<7 


An 


HZ 


1 22 

^n 5 a 


487T 

-\ T a T bc + T^TT ab + - (4) 


12 


_ _ cd _ 

9abT T c d 


_ 3 _- 7 - 2 ( 4 ) 

24 




9 ab ? 
e (5) T 




SD-r, a I , b 


(5) T BD n B n D + -e (5) T 


(4) 


9 ab 


( 8 ) 

( 9 ) 


( 10 ) 


( 11 ) 


Here we stress that the effective 4-dim cosmological constant can be set to 
zero in the present case of an extra timelike dimension by properly fixing the 
bulk cosmological constant as A 5 = a 2 . It is important to notice that 

for a 4-dim brane embedded in a conformally flat empty bulk we have the 
absence of the Weyl conformal tensor projection E ab , and of F ab in Eq. (J7J). 

On the other hand, Codazzi’s equations imply that 

V a K - V„At = ie/c|v 4 rj . (12) 


By imposing that V&t£ = 0, the Codazzi conditions read 


W a E ab = 4\7 a n ab + V a F ab . (13) 

where V a is the covariant derivative with respect to the induced metric ^g ab - 
Eqs. © and (1T3j) are the dynamical equations of the gravitational field on 
the brane. In the following section we drop the index (4) in the geometrical 
quantities on the brane. 


3 The model 

Let us consider a Bianchi IX spatially homogeneous geometry on the four- 
dimensional brane embedded in a five-dimensional, conformally flat and empty 





bulk (E a b = 0 = F a b) with a timelike extra dimension (e = 1). In comoving 
coordinates on the brane, the line element can be expressed as 

ds 2 = dt 2 - ( 9 l ) 2 - ( 9 2 ) 2 - (0 3 ) 2 , (14) 


where t is the cosmological time and 

9 1 = M(t) u 1 , 9 2 = N(t ) (u 2 , 6 3 = R(t) cu 3 . 


(15) 


Here M(t), N(t) and R(t) are the scale factors of the model and the uj 1 
{i = 1,2,3) are Bianchi-type IX 1-forms satisfying 

dJ = ^e ijk u j A oo k , (16) 

where d denotes the exterior derivative. 

The matter content of the model is assumed to be dust, whose energy 
density p is measured by the comoving observers with 4-velocity u a = 5q. By 
imposing that the energy-momentum tensor of dust, 

r ab = pu a u b , (17) 


is conserved separately, namely V a r ab = 0, we obtain 

Co 

P MNR 1 


(18) 


where Co is a constant of motion connected to the dust energy. The compo¬ 
nents of tensor n a ^ are given by 


Ilab ^2 P 9abi 


(19) 


so that Codazzi’s equations, 


v fl n ab = 0, (20) 

are identically satisfied. Therefore, Eq. (J7J) reduces to 

C a b T A (jab 8llGj\jT a b Kc^Tlabi (^l) 

which are the modified Einstein’s field equations for the model. As the 
Gauss-Codazzi equations are automatically satisfied via (12UD and (1211) . the 
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assumption of a conformally flat empty bulk is consistent. We also see that 
as e = 1 (a timelike extra dimension) the term II a b in (121|) acts as a potential 
barrier to the dynamics avoiding the singularity. 

In terms of the metric functions (jT4|) equations (1211) correspond to the 
modified Friedmann’s equations of the model, having a first integral that can 
be expressed as the Hamiltonian constraint 


H 


1 ( M 2 N 2 R 2 2 2 2 \ 

8 V Wr Pm ~ ~MR Pn ~ ~MN Pr + m PnPr + n PmPr + r P mPn ) 

---[M 4 + N 4 + R 4 - (M 2 - N 2 ) 2 - (R 2 - M 2 ) 2 - (R 2 - IV 2 ) 2 ! 

2MNR V ; V ; V J 1 

E 2 

2AMNR - 2E 0 + k 2 = 0 (22) 

MNR v ’ 


where pm , Pn and p R are the momenta canonically conjugate to M, N and 
R, respectively. E 0 = 8-kGnCq and k 2 = (SttGat) - 1 ^!" 1 . From Hamilton’s 
equations we obtain the following dynamical system 


M 

N 

R 


Pm 


Pn 


d H 1 (p N i p R M __ ^ 

dp M ~4:\R + N NR PM h 
dH _ 1 / p M p R N \ 

dp N ~ 4 V R + M MR PN ) ’ 
dH _ 1 /p M Pn C \ 

dp R ~ 4 V N + M MN Pr ) ’ 
dH _ 1 (p 2 M N 2 R 2 | 2 \ 

dM 8 \NR M 2 R Pn M 2 N Pr+ M 2PnPr ) 

+ 2 AHNR ^ + N4 + R4 ~ ( R2 ~ N2 ) 2 “ “ M2 ) 2 _ ( M2 “ iy2 ) 2 ] 

+2ANR ~ mm lM3 + M(R2 ~ m2) ~ M{M2 _ ^ 

dH l ( P% iW 2 i? 2 , 2 \ 

= 8 VAFR “ ~~ MN 2Pr + Jp P mPr ) 

+ 2MN 2 R lMi + N * + R4 ~ _ Ar2 ) 2 _ “ M2 ) 2 - ( M2 - iy2 ) 2 ] 

+2AVB - ^[JV 3 + iV(i? - JV 2 ) - ]V(JV 2 - M 2 )] + 
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dH If pi M 2 Nr, 2 \ 

Pn “ ~~M ~~ 8 VMiV ~~ M Pm ~~ MB? Pn + & PMPN ) 

+ —4tb?[M 4 + AT 4 + i? 4 - (i? 2 - iV 2 ) 2 - ( R 2 - M 2 ) 2 - (M 2 - iV 2 ) 2 ] 
1MJM H z 

+2AMN ~ mm lR3+R{M2 ~ R2) ~ R(R2 ~ Af2)l + mm - (23) 

Equations (|22|) and (123|) are equivalent to the modified field equations f[2T]) . 


4 The structure of the phase space 


In this section we will examine the basic structures that organize the dynam¬ 
ics of the system in the phase space. The first of these are the set of critical 
points of the system given, from Eqs. ([23]) ). by M = N = R = Mo and 
Pm = Pn = Pr = 0, where M 0 satisfies the equation 


Mn 


if + r2 

4A 


4 = 

2A 


(24) 


We can observe that the critical points, determined by the positive real roots 
of (El) , depend on their respective critical energy appearing in the third 
term of the left-hand-side of the equation, as a consequence of the bulk-brane 
interaction. 

We must also consider the further relation 

l M ° +k2 l§" 2AM ° 3 - 2Ecr= °> ( 25) 

obtained by evaluating the Hamiltonian constraint (|22l) at the critical points. 
Solving ([25D f° r Ecr we will restrict ourselves to the root 


E 


cr 


Ml 

K 2 




3k 2 
2 M 2 


+ 2n 2 A ) 


(26) 


which yields the correct result in the general relativity limit [39] (ft 2 —> 0 or 
equivalently |cr| —>■ oo). Combining Eqs. (|24l) and (]26l) we obtain for the 
critical points the two real positive solutions 


M 


0l,2 


K (3 ± y/1 - 16k 2 A ) 

2(1 + 2 k 2 A) yj(l- 4k 2 A ± Vl - 16K2A )/(l + 2 k 2 ~A) 


(27) 
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with M 0l < M 02 . The equality occurs for A = 1/16k 2 , the case of just one 
critical point; for A > 1/16k 2 no critical points exist. In the following we 
are going to restrict ourselves to the case A < 1/16k 2 . As we will see, this 
condition is necessary for the presence of homoclinic orbits that establish the 
chaotic behavior of the dynamics. The respective energies associated with 
the critical points are obtained by substituting M 0l 2 in yielding 

r _ « ( 3 ± \/l - 16 k 2 A ) 2 , 00 ^ 

^Cr/i 21 /- . 

8 (1 + 2k 2 A) 3 / 2 \/1 — 4k 2 A ± \/l — 16 k 2 A 

Much of our understanding of nonlinear systems derives from the lin¬ 
earization about critical points and from the determination of existing in¬ 
variant submanifolds, which are structures that organize the dynamics in 
phase space. The system under examination here presents a two-dimensional 
invariant manifold of the dynamics defined by 

Pm = Pn = Pr , M = N = R. (29) 

This invariant plane is actually the intersection of two four-dimensional in¬ 
variant submanifolds, defined by (M = N, pm = Pn) and (N = R, pn = Pr)- 
The critical points obviously belong to the invariant plane. 

Finally a straightforward analysis of the infinity of the phase space shows 
that it has two critical points in this region, one acting as an attractor (stable 
de Sitter configuration) and the other as a repellcr (unstable de Sitter con¬ 
figuration). The scale factors M, N and R approach the de Sitter attractor 
as M = N = R ~ exp(^/A/3 t), so that the two de Sitter configurations 
also belong to the invariant plane. The phase picture of the invariant plane 
is displayed in Fig. [U in the variables (x,p x ) defined in section [5] 

To proceed let us now linearize the dynamical equations (125|) about the 
critical points (M = N = R = M 0 ., p M = p N = p R = 0), i = 1, 2. Defining 

X = (M-MqJ, W = (p M - 0), 

Y = (N — M 0i ), K=( Pn - 0), (30) 

Z = (-R —M 0i ), L = (p R - 0), 
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we obtain 


where 


/ 


V 


X \ 


/ 0 0 0 —a a a \ 


( X \ 

Y 


0 0 0 a —a a 


Y 

Z 


OOOcx a —a 


Z 

w 


/3 7 7 0 0 0 


W 

K 


O 

O 

O 

C“~~ 


K 

L ^ 


0 

0 

0 


\ L ) 


a = 


0 = 


4 M 0i 

7 = 2 A M 0i - 


M 0i 

3 

2Mo ; 


2 « 2 ^ (0 

^ : 


k 2 u 2 


r W 




The associated characteristic polynomial results 

P( A) = (A — + (3a )(A + \/2^ol + /3a ) 

x(A - ^2a(j- P) ) 2 (A + y 2 a (7 - P) ) 2 , 


( 31 ) 


(32) 


(33) 


with roots 


x (i) 

Ap) 


±i 


.V2 


M< 


Oi 


± 



1 

2 ’ 


(34) 

(35) 


where (l24|i was used. The pair of imaginary eigenvalues (|34|) has multiplicity 
two, characterizing a center-center structure. The analysis of the center- 
center structure will reveal a manifold of linearized unstable periodic orbits 
with the topology of S' 3 . The extension of this manifold to the nonlinear 
domain constitutes the center manifold [211 25| °f unstable periodic orbits, 
parametrized with the constant of motion Eq (with Eq < E CI ), which will 
play a central role in our discussions in the next section. 

Using (l25|i and (1271) one can show that the pair of eigenvalues (|35|) are 
imaginary for the critical point i = 1 and real for the critical point i = 2 . As 
we shall see below, we have (in the latter case i = 2 ) that the critical point P 2 
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is a saddle-center-center about which the dynamics has the topology R x S 3 . 
On the other hand, the critical point Pi is a center-center-center critical 
point, about which the dynamics has the topology S 1 x S 3 corresponding to 
perpetually oscillatory Bianchi IX universes. 

Finally we should note that, in the limit case of a single critical point 
(when 16 k 2 A = 1), the second pair of eigenvalues (13511 are zero and no saddle 
structure is present in the dynamics. The analysis of this case will not be 
undertaken here. In the remaining of this section our discussion follows the 
lines of [39] done for the general relativity case. 

To display the structure of the linearized motion, we start by diagonaliz¬ 
ing the linearization matrix of (13T]) with the use of a similarity transformation 
whose columns are composed of six independent eigenvectors of the lin¬ 
earization [ID]. A judicious choice of 9ft yields primed variables defined by 
the transformation 

x' = i(x + r + z), 

y' = — (x-y), 

M 0i K h 

z' = -^-(X + Y-2Z), (36) 

M o i 

W' = (W + K + L ), 

A" = ^(W-K), 

L' = (w + K-2L). 

6 

In these new variables, the quadratic Hamiltonian about the i-th critical 
point is expressed in the form 

„ = I(i^-6 9 ^)-(P J x' 2 + M 0 ,y' 2 ) 

- (2k L ' 2 + ^ Z ' 2 ) +2{E ^~ E ° ] ' ^ 

where 


q = 6(A M 0i - 


K 


■E 


JO 


r (i) 


(38) 
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These primed variables are conjugated to the pairs according to [A"', W'] = 
1, \Y',K'] = 1, [Z',L'] = 1, other Poisson brackets (PB) zero. The Hamil¬ 
tonian (1371) is separable, and we can identify the following constants of the 
linearized motion 

E q 

= j(jV- 6 ,V 2 ), 

(39) 

E r otl 

“ 2 M 0 3 / V ' ’ 

(40) 

E r ot2 

= 3 . L' 2 + M °‘ Z' 2 , 

2 Ml 3 

(41) 

Qi 

= AIoi Y'Z' + 1 , K'L', 

3 2 M 3 0i 

(42) 

Q 2 

1 (l'Y’ 1 K'Z '^], 

2 M 0i V 3 P 

(43) 


in the sense that they all have zero PB (137|) . The first three constants appear 
as separable pieces in the Hamiltonian fj37)l . 

The case of E q demands a separate analysis for the two critical points. 
From previous relations we have that q > 0 for the critical point P 2 , so 
that E q corresponds to the energy associated with the motion in the saddle 
sector. We remind that this is connected to the fact that the second pair of 
eigenvalues (1351) are real for P 2 . For the critical point P± in which q < 0, E q 
corresponds to the rotational energy in the additional rotational sector of the 
dynamics about Pi which has the structure S 1 x S 3 as mentioned already. 
The center-center-center critical point Pi corresponds to a stable Einstein 
universe configuration with perpetually oscillatory orbits in its neighborhood. 

In the following our focus will be the dynamical phenomena connected to 
the presence of the saddle-center-center critical point P 2 in the phase space 
of the model. We remark however that the analysis of the center manifold of 
unstable periodic orbits can also be applied to the case of the critical point 

Pi- 
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5 The dynamics about the saddle-center-center 
critical point 


We will now proceed to describe the topology of the general dynamics in the 
linear neighborhood of the saddle-center-center critical point P 2 for which 
q > 0 . 

If E q = 0 two possibilities arise. The first possibility is W' = 0 = X'. 
The total energy in this case is E rot 1 + E rot2 , the sum of the energies of the 
rotational motion in the linear neighborhood of the center-center manifold, 
corresponding to the motion on 2-dim tori [21] . The remaining two constants 
Q 1 and Q 2 are additional symmetries that arise due to the multiplicity two of 
the imaginary eigenvalues and are connected to the fact that the linearized 
dynamics in the center-center sector is that of a 2 -dim isotropic harmonic 
oscillator. They are not all independent but related by 

4 E ro t 1 E r0 f2 12 Q\ + 6 Ql (44) 


The motion in the constant energy surfaces W' = 0 = X' are periodic orbits 
of the 2-dim isotropic harmonic oscillator, with Hamiltonian 



— 2(i? cr(i) — E 0 ) — 0. 


(45) 


The above equation shows that Eq — E cr < 0 is necessary for the dynamics 
in the rotational sector, defining a condition for the existence of the center- 
center manifold of periodic orbits. 

By a proper canonical rescaling of the variables in (ITjjl we can see that 
these constant energy surfaces are hyperspheres and that the constants of 
motion Q 1 , Q 2 and Q 3 = E rotl —E rot2 satisfy the algebra of the 3-dim rotation 
group under the PB operation, namely, 


[Qi, Qj] — ^ k Qk- 


(46) 


The constant of motion Q 1 considered as a generator of infinitesimal contact 
transformations has a peculiar significance in characterizing the topology of 
the underlying group of the algebra (l46lh While Q 2 generates infinitesimal 
rotation of the orbits, Q\ generates infinitesimal changes in eccentricity. The 
action of Qi is to take an orbit - let us say nearly circular - and to transform 
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it into an orbit of higher and higher eccentricity until it collapses into a 
straight line. Continued application of Qi produces again an elliptic orbit, 
but now traversed in the opposite sense, so that it takes a 720° to bring the 
orbit back into itself. The two-valuedness of the mapping arises from the 
fact that the orbits are oriented. Therefore the group generated by these 
constants of motion is homomorphic to the unitary unimodular group [TT] so 
that the topology of the center-center manifold is in fact S 3 . 

Due to the separate conservation of E rotl and E rot2 (cf. (H5jl ) one can show 
that the center manifold in the linear neighborhood of the critical points is 
foliated by Clifford 2 -dim surfaces in S 3 [42], namely, 2 -tori 5s e 0 contained 
in the energy surface E 0 = const. Such surfaces, as well as the S 3 manifold 
containing them, depend continuously on the parameter E 0 . We remark that 
these two tori will have limiting configurations E rotl = 0 or E rot2 = 0, and 
correspond to the case of maximum eccentricity (for instance, a straight line 
in the plane (Y r , Z')). 

The second possibility to be considered is W' = ±6y fqX' . It defines 
the linear stable Vs and unstable Vjj manifolds of the saddle sector. Vs 
and Vjj limit regions I (E q = Eh yp < 0) and regions II (E q = Eh yp > 0) 
of motion on hyperbolae which are solutions in the separable saddle sector 
Eh yp = W' 2 — 6qX' 2 ). Note that the saddle sector depicts the structure of 

the neighborhood of P 2 in Fig. 1 , with Vjj and Vs tangent to the separatrices 
at P 2 . The direct product of 5s Eo with Vs and Vjj generates, in the linear 
neighborhood of the critical point (i = 2 ) the structure of stable (5s Eo x Vs) 
and unstable (5s e 0 x W/) 3-dim tubes which coalesce, with an oscillatory 
approach to the tori 5s e 0 for t —y oo. The energy of any orbit on these 
tubes is the same as that of the orbits on the tori 5s e 0 - These structures are 
contained in the 4-dim energy surface H = E 0 such that (Eq — E cr ) < 0. 
We should recall that the tubes constitute a boundary for the general flow 
and are defined by E q = 0 in the linear neighborhood of the critical point. 
Depending on the sign of E q the motion will be confined inside the 4-dim tube 
(for E q < 0) and will correspond to a flow separated from the one outside 
the tube (for E q > 0). A detailed examination of the above motion and its 
extension to the nonlinear regime will be done in the next section. 

The extension of our analysis beyond a linear neighborhood of critical 
points could be made by implementing normal forms [43l [44] and associated 
coordinates, modulo their radius of convergence. We will instead propose 
here a suitable canonical transformation which will allows us to obtain an 
exact analytical form for the center manifold as well as a sufficiently accurate 
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description of the phase space dynamics in extended regions away from the 
critical points. In particular we can examine the behavior of the nonlinear 
extensions Ws and Wjj of, respectively, the linear stable (3\g 0 x Vs) and linear 
unstable (Xse 0 x Vu) manifolds [29] emanating from the neighborhood of the 
saddle-center-center P 2 . Let us introduce the canonical transformation with 
the generating function 


G = ( MNR) 1/3 p x + ^ p y + —j^rPz, 
where p x , p y and p z are the new momenta, resulting in 


and 


x = ( MNR ) 113 , y = 


M 

IV’ 


z = 


MN 


Pm 

Pn 

Pr 


1 NR 1 ^ 
3 ( MNR) 2 /z Px + N Py 


N 


MR 


3 (MNR) 2 / 3 
1 MN 
3 (MNR) 2 / 3 


M 

N 2 ' 1 


R 2 

M 


Pz , 


P x at2Pv + N 2 ^ Z1 


Px ~ 


2 MN 


R 3 


-Pz- 


(47) 


(48) 


(49) 


Here, the variable x is obviously the average scale factor of the model. In 
these new canonical variables the equations of the invariant plane reduce to 


V = 1, z = l, p y = 0 = p z . (50) 

It is then clear that (x,p x ) are variables defined on the invariant plane. In 
these variables the phase space picture of the invariant plane is given in Fig. 
CD The separatrices S emerging from the saddle-center-center P2 separate the 
invariant plane in three disconnected regions, region / of oscillatory universes 
and regions II and III of one bounce universes. They are constituted of 
three branches, namely, the separatrix that divides the regions / and II and 
makes a homoclinic connection with the critical point P 2 and two others 
that approach the de Sitter asymptotic configurations for t —> ±00. The 
first branch will play a fundamental role in our following discussions and 
will be referred to as separatrix, except where a qualification is needed to 
avoid confusion. The center-center-center P\ corresponds to a stable Einstein 
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Figure 1 : The phase portrait of the invariant plane for k 2 = 0.5/Eq and 
A = 0.01. The critical points Pi and P 2 belong to the invariant plane. The 
separatrices S emerging from P 2 separate the plane in three disconnected re¬ 
gions, one of oscillatory universes the other two of one bounce universes. The 
separatrix dividing the regions / and II makes a homoclinic connection with 
P 2 . The linear neighborhood of P 2 depicts the motion of the saddle sector 
with Vs and V\j tangent to the separatrix at P 2 . The figure is constructed in 
the canonical coordinates defined in (l47|h 
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universe configuration that occurs due to the bulk-brane interaction term 
proportional to Eq in the Hamiltonian (I22|) . 

In the variables (x,p x ,y,p y , z,p z ) the full Hamiltonian (l22]i assumes the 
form 


1 


2 y 2 2 


H(x, PxiViPy, Z, p z \ E 0 ) = — Px - —p y 
3 P 


+ 


2x 3 

xy 

r l/3 


Pi 


X 


2z 4 / 3 


-xz 2 ' 3 y 2 
2 y 


1 

V 


■ xz 2/3 + 


X 


yz 


1/3 


K 2 F 2 

+ xz 2/3 + 

X 6 


2 Ax 3 - 2A n = 0. 


(51) 


We remark that the linearization of fj48ll - fj5T]l about both critical points 
Pi and P 2 of the dynamical system (1231) yields exactly the transformation 
(J36]), and that the variables ( y,p y ,z,p z ) correspond to the primed variables 
(A'', Y' : L ', Z') defined on the center-center manifold S 3 about a linear neigh¬ 
borhood of P 2 - 

The new canonical variables are most convenient since they separate the 
degrees of freedom of the system into pure rotational modes, ( y,p y ) and 
(z,p z ), and the expansion/contraction mode {x,p x ) connected to the invari¬ 
ant plane. This can be illustrated by implementing the expansion of the 
dynamical system generated from (l5Tf about a linear neighborhood of the 
invariant plane, producing a linearized Hamiltonian parametrized by the vari¬ 
ables (x(t),p x (t)) describing the curves in the invariant plane, for instance, 
in the region / of periodic orbits bounded by the separatrix S homoclinic to 
P 2 . This is analogous to the usual expansion of a dynamical system about a 
periodic orbit. Using (15TJ) . we then obtain 


where 


H 


£inv ~ 2 ^^ + 3p ^ ~ X<kV ~ 1 ^ 
^x(z - l) 2 = 2 E 0 , 


£u 


1 2 3x k 2 Aq 

2Ax Px 2 x 3 


— 2 Ax 3 = const. 


(52) 


(53) 
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The resulting dynamical equations are 


8y = 

3 fiPyi 

x 6 

Py = 

2x 8y , 

5z = 

3 A 

- 3 Sp z , 

X 6 

8p z = 

2x 

y 0Pv, 


( 54 ) 


where 5y = (y - 1), Sz = (z - 1), 8p y = (p y - 0), and 8p z = (p z - 0). 
The linearization matrix of (|54l) has imaginary eigenvalues A = ±i\/2/x(t), 
both with multiplicity two, corresponding to elliptical modes in the linear 
neighborhood of the invariant plane so that the motion is oscillatory about 
the invariant plane. 


6 The non linear center manifold and the ho¬ 
moclinic cylinders 

The nonlinear extension of the center manifold, by continuity, maintains the 
topology S 3 but it can no longer be decomposable into E rotl and E rot2 so 
that now only the 4-dim tubes with the topology R x S 3 are meaningful 
for the nonlinear dynamics. Similarly the extension of the structure of the 
4-dim tubes away from the neighborhood of the center manifold are to be 
examined, and our basic interest will reside in the stable and unstable tubes, 
Ws = Vs x S 3 and Wjj = V v x S' 3 , that leave this neighborhood. The 
tubes have the structure of 4-dim spherical cylinders (of co-dimension 2), 
one less dimension than the energy surface, and act therefore as separatrices, 
separating the energy surface in two dynamically disconnected parts. The 
2-dim invariant plane, defined by (|29l) . is contained in a 6-dim phase space 
and it is obvious that, contrary to examples in lower dimensional systems, it 
does not separate the dynamics in disjoint parts. In fact the general motion 
about the curves of the invariant plane is an oscillatory flow confined in the 
interior or exterior of 4-dim tubes R x S 3 , so that the invariant plane (or 
more properly, one of the curves of the invariant plane) can be thought as a 
structure in the center of the tubes. This latter fact is of crucial importance 
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Figure 2: The sections (y = 1 ,p y = 10~ 3 ) (left panel) and (z = 1 ,p z = 
0) (right panel) of the center manifold for five values of E 0 . We note the 
deformation of the center manifold S 3 in the nonlinear domain as (E cr — E 0 i ) 
increases. Here E cr = 2.5127254138199464. 


in the discussion of the transversal crossing of the 4-dim cylinders Ws and 
Wjj made in section [0 

The nonlinear extension of the center manifold in the canonical variables 
(y, z,py,p z ) is obtained by substituting x = x cr and p x — 0 in (l5T]h yielding 
after some manipulation the exact analytical expression 


H r = 


?>z 2 


2k P * + ^k P * + Xcr h + 2^ 


-z 2 'V + Z 


cr 

2/3 


y 


2 y 2 


yz 


1/3 


-1/3 


,2/3 


,E n 2 


K 


x-i 


— 2 (E cr — E 0 ) — 0, 


(55) 


where x cr and E cr are respectively the coordinate and the energy of the crit¬ 
ical point /" 2 • The form (1551) adopted above for the center manifold equation 
makes explicit its dependence on the parameter (E cr — Eo). For E$ = E cr the 
center manifold reduces to the critical point. The domain of E$ defining the 
center manifold satisfies the constraint to E 0 < E cr , as already discussed; the 
case of the linear version (147>|) corresponds to ( E cr — E 0 ) sufficiently small. 
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Figure 3: (left panel) The wireframe shows the 2-dim section p y = 10 “ 3 of 
the center manifold in the space ( y,z,p z ) for Eq = 2.512725. The solid line 
indicates an orbit for 0 < t < 2000 with initial conditions (xO = xcr, r/0 = 
1,^0 = l,px0 = 0,py0 = 10” 3 ,px0 = 0.008177970306403768) on the center 
manifold, (middle panel) An expanded piece of the center manifold on the 
left, including the same orbit, (right panel) The wireframe shows the 2- 
dim section p y = 10~ 3 of the center manifold in the nonlinear domain, with 
Eq = 2.3. The solid line corresponds to a one-bounce orbit which moves 
towards large values of z (when t ~ 37, z(t) ~ 63) before escaping to the de 
Sitter attractor when t — 70. The initial condition for this orbit, (xo = x cr , 
y 0 = 1, z 0 = 1, p xo = 0, p yo = 10^ 3 , p ZQ = 5.879276478889785) is taken on 
center manifold. This orbit remains on the center manifold for a time up to 
t ~ 12. Here E cr = 2.5127254138199464. 

As (E cr — E 0 ) increases we have a nonlinear center manifold parametrized 
by the energy Eq. In general the center manifold is a 3-dim submanifold of 
the 6-dim phase space contained in the 5-dim energy hypersurface H = Eq. 
In Figs. | 2 ] we plot the sections (y = 1 ,p y = 10~ 3 ) and (z = 1 ,p~ = 0) of 
the S 3 center manifold (155|) showing its deformation in the nonlinear regime 
as the values of ( E cr — Eq) increase. We adopted the values A = 0.01 and 
k 2 = 0.5 so that the associated critical energy E cr = 2.5127254138199464 
and x cr = 4.974148895632555 for the saddle-center-center P 2 . In the Figures 
we selected five values for Eq. 

As we have already seen the canonical coordinates {y,p y ,z,p z ) cover the 
center manifold S 3 and therefore we will use them not only to examine the 
stability of the motion restricted to S 3 but also to obtain an accurate descrip- 
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Figure 4: (left panel) The wireframe shows a piece of the 2-dim section 
p z = 0 of the center manifold in the coordinates (x, y,p y ) for E 0 = 7.9096539. 
ffere we fixed the parameters A = 0.001 and n 2 = 0.5, with corresponding 
E cr it = 7.909653935312942. The solid line is an orbit of the full dynamics for 
0 < t < 645 with initial condition (x 0 = x cr(2) , y 0 — 1, z 0 — 1, p XQ = 0, p yo = 
0.023600434934712197, p ZQ = 0) taken on the center manifold. The orbit 
remains on the center manifold up to t ~ 500 when it leaves its neighborhood 
and is driven towards its first bounce, (middle panel) The same orbit of the 
left figure for 652 < t < 1320. We can see that the orbit returns from 
its first bounce to a sufficiently small neighborhood of the center manifold 
for t ~ 652, before escaping to its second bounce when t ~ 1320. (right 
panel) The solid line indicates the same orbit of the previous figures but for 
1330 < t < 2020. This orbit returns from its second bounce to a sufficiently 
small neighborhood of the center manifold for t ~ 1330 before escaping to 
the de Sitter attractor when t ~ 2020. 
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Figure 5: Plot of p y (t) and y(t) corresponding to the orbit of Figs. ® showing 
the typical oscillatory behavior of the rotational modes in the dynamics of the 
model. We note the decrease of the amplitude of p y (t) and the amplification 
of the conjugated y{t) in a neighborhood of the bounce. This pattern is 
analogous for the other variables of the orbit (p z , z ). 

tion of whole the dynamics (x(t),p x (t),y(t),p y (t),z(t),p z (t)) emerging from 
a neighborhood of the center manifold. In the following we will numerically 
illustrate this behavior. We must remark that we do not make use here of 
the displacing (in the direction of the unstable cylinder) of initial conditions 
taken on the invariant center manifold, as the shooting method in [34], but 
instead we make use of the instability of the motion on the center mani¬ 
fold which computationally conserves the Hamiltonian constraint (T5TT) for all 
t. Actually in all our numerical simulations the error in the Hamiltonian 
constraint (1511) is checked to remain < 10~ 13 for the whole computational 
domain. 

To start let us fix the parameters n 2 = 0.5 , A = 0.01 as in Figs. [2] In 
Figs. [3] we now show the 2-dim sections p y = 10~ 3 of the center manifold 
for Eq = 2.512725 (left). The solid line indicates an orbit with initial condi¬ 
tions obviously satisfying (15511 . This orbit is evolved with the full dynamics 
generated from the Hamiltonian (]5TT) and remains on the center manifold for 
0 < t < 2000. A piece of this center manifold is displayed in Fig. [3] (middle) 
where the solid line describes the same previous orbit. Fig. [3] (right) displays 
the section p y = 10~ 3 of the center manifold for E 0 = 2.3. The solid line cor¬ 
responds to a one-bounce orbit which moves towards large values of z (when 
t ~ 37, z(t) ~ 63) before escaping to the de Sitter attractor when t ~ 70. 
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The initial condition for this orbit is taken on the center manifold S' 3 ; it 
remains on the center manifold up to £ ~ 12. This increase of the dynamical 
instability is actually due to the large value of (E crit — E 0 ) ~ 0.213, causing 
the orbit to leave rapidly the center manifold towards the bounce, satisfying 
however the exact dynamics within an error < 10~ 13 . 

A second set of experiments is displayed in Figs. 0 where we examined 
the oscillatory motion originating in the 2-dim section p z = 0 of the center 
manifold for A = 0.001, k 2 = 0.5 and E 0 = 7.9096539. The continuous solid 
line shown in the left panel is an orbit of the full dynamics for 0 < t < 645 
with the initial condition taken on the center manifold. This oscillatory orbit, 
which is initially periodic, remains on the center manifold up to t ~ 500 when 
it leaves this neighborhood and is driven towards its first boimce. In the 
middle panel the same orbit is shown for 652 < t < 1320. We see that the 
orbit returns from the first bounce to a sufficiently small neighborhood of the 
center manifold for t ~ 652 before escaping to its second bounce at t ~ 1320. 
In the right panel we have the same orbit of the previous figures but for for 
1330 < t < 2020. This orbit returns from its second bounce to a sufficiently 
small neighborhood of the center manifold for t ~ 1330, before escaping to 
the de Sitter attractor when t ~ 2020. These numerical simulations also 
reveal a typical behavior of the dynamics as we decrease Eq. In fact the 
increase of ( E crit — E 0 ) makes an orbit, with initial conditions taken on the 
center manifold, to rapidly leave this neighborhood indicating a dynamical 
instability (albeit the accuracy of the exact dynamics) as shown in Fig. [3] 
(right). 

We remark that the oscillatory behavior of the orbit in the phase space 
sectors ( y,p y ) and ( z,p z ) is typical, even when the orbit tends asymptotically 
to one of the deSitter attractors. This is illustrated in Figs. 0 where we plot 
the time behavior of p y and y of the orbit discussed in Figs. [H We note a 
decrease of the amplitude of p y and an amplification of the amplitude for the 
conjugated y in a neighborhood of the bounce, this pattern is analogous for 
the other variables ( p z , z) of the orbit. 

Finally we give a numerical illustration of the stable and unstable cylin¬ 
ders emanating from the center manifold which are a nonlinear extension of 
the ^seq'xVs and Ag 0 x Vu, with Ae 0 C S 3 defined in a linear neighborhood 
of the saddle-center-center P 2 . We must recall that these cylinders are actu¬ 
ally composed of orbits that have the same energy ( E cr — Eq) of the center 
manifold and coalesce to it as t ± 00 . In Fig. [6] we display the stable Ws 
and unstable Wu cylinders emanating from the neighborhood of the center 
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Figure 6: A numerical illustration of the unstable cylinder Wjj (black) 
spanned by 16 orbits, with initial conditions taken on a circle in the do¬ 
main ( y,p y ) of the center manifold for E 0 = 7.9096 , emerging towards 
the bounce. From the same domain of initial conditions the stable cylin¬ 
der Ws (gray) emerges towards the bounce. The projection of the figure 
in the plane ( x,p x ) “shadows” the separatrix of the invariant plane. Here 
E crit = 7.909653935312942. 

manifold towards the bounce, guided by the separatrix dividing the regions 
I and II of the invariant plane. We emphasize that the separatix guiding the 
cylinders is actually a structure inside the cylinders. We fixed the parame¬ 
ters A = 0.001 and n 2 = 0.5, as in Figs, [dj and took Eq = 7.9096 so that 
(E cr - Eq) ~ 10 -5 . 

A comment is in order now. Since the cylinders W$ and Wjj are 4- 
dim surfaces they obviously separate the 5-dim energy surface defined by 
the Hamiltonian constraint (|5TT) in two dynamically disconnected pieces, a 
fact that will be fundamental in the characterization of chaos in the case 
of an eventual transversal crossing of Ws and HA[291 [22] . We remark that 
in Fig. [H] the projection of the cylinders on the plane ( x,p x ) “shadows” 
the separatrix of the invariant plane, as expected since the separatrix is a 
structure contained in the interior of the two tubes. As the separatrix in 
question makes an homoclinic connection to the saddle-center-center critical 
point P 2 this fact leads necessarily to the transversal crossings of Ws and 
Wu, a dynamical phenomenon that we examine in the next section. 
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7 The transversal crossings of the cylinders 
and the homoclinic intersection manifold: 
a chaotic saddle 

The results of the previous sections showed that two 4-dim cylinders, one sta¬ 
ble Ws = RxS 3 and one unstable Wu = RxS 3 , emerge from a neighborhood 
of the saddle-center-center P%. The center manifold S 3 is the locus of the ro¬ 
tational degrees of freedom of the phase space dynamics and is parametrized 
with the energy E 0 (E 0 < E crit . ft encloses the critical point P 2 and tends 
to it as Eq — > E crit . At this limit the cylinders Ws and Wu reduce to the 
separatrix S, which makes an homoclinic connection of P 2 to itself in the 
invariant plane. The separatrix is a structure inside the cylinders, about 
which the flow with the oscillatory degrees of freedom (y,p y , z,p z ) proceeds, 
guiding the cylinders towards the bounce (cf. Fig JT]) and leading to their 
eventual crossing. The first crossing is expected to occur in a neighborhood 
of the bounce {x = Xb,p x = 0), where Xb is the scale factor of the bounce for 
the orbits at p x = 0. In order to detect this first intersection we will adopt 
as the surface of section [45] the 4-dim surface E : (x — Xb,p x = 0). This first 
transversal crossing of the cylinders will be the main object of the present 
section. 

Due to the conservation of the Hamiltonian constraint (1411) we have that 
at the bounce 


H(x = x b , p x = 0, y, p y , z, p z ; Eq) = 0, (56) 

which is the equation of a closed surface with the topology of S 3 . The 
transversal crossing of the stable Ws and unstable Wjj cylinders at the bounce 
will therefore be a set of points contained in the transversal intersection of 
two 3-spheres defined by (144]) . then a S 2 . These points define a set of orbits 
that are contained both in the stable and the unstable cylinders and are bi- 
asymptotic (homoclinic) to the center manifold S 3 (Eq). They are denoted 
as homoclinic points and homoclinic orbits. Therefore the set of homoclinic 
points has the compact support S 2 . The presence of a homoclinic orbit in 
the dynamics is an invariant signature of chaos in the mo del [271 13 Ij . The 
homoclinic intersection manifold has the topology RxS 2 and consists of 
all homoclinic orbits biasymptotic to the center manifold. In this sense, a 
chaotic saddle\ 34] associated with S 3 is defined, indicating that the set of 


intersection points of the cylinders has the nature of a Cantor-type set with 
a compact support S' 2 [515]. 

A complete numerical study of the intersection of the 4-dim cylinders Ws 
and Wu is beyond the scope of the present paper (it will be considered as 
the subject for a future publication). Here our numerical experiments will 
be restricted to the dynamics on the two 4-dim invariant submanifolds of 
the 6-dim phase space defined by (i) M = N, p M = p N (or equivalently 
y — 1, p y — 0), and (ii) N = R, p^ = pr (or equivalently y = z, p y = 3 p z ). 
The denomination invariant submanifolds derives from the fact that each 
of them is mapped on itself by the general Hamiltonian flow (l23]h in other 
words, invariant under the flow. We will then examine the intersection of 
2-dim stable and unstable cylinders in these two 4-dim invariant manifolds, 
according to the following experiments. 

To start we fix the parameters n 2 = 0.5, A = 0.001, with corresponding 
E cr = 7.9096539353149939 and x cr = 15.8034567969528. The total energy of 
the system is taken E 0 = 7.9096, so that the energy available to the rotational 
modes will be given by ( E cr — E 0 ) ~ 10 -5 . 

In the first experiment we take (xo = x cr n, p xo = 0), and fix initial 
conditions on the 4-dim invariant submanifod (i), namely, with (yo = 1, 
Pyo = 0)- That is, we are restricting ourselves to a particular domain of 
initial conditions in the sector ( z,p z ) of the center manifold S 3 , which has 
the topology of S 1 and is defined by the Hamiltonian constraint (El]) as 

H(x = x cri p x = 0,y = l,p y = 0 ,z,p z ) = 0. 

By performing the evolution of orbits from a large set (of the order of 1, 500) 
of initial conditions in the above domain, the exact dynamics actually evolves 
a 4-dim invariant subset (x, p x , z, p z ) of the full 6-dim phase space as expected 
due to our restriction to the 4-dim invariant manifold (y = l,p y = 0). In this 
particular experiment, we have that under the exact dynamics no motion is 
present in the sector (y,p y ). We generate one 2-dim stable and one 2-dim 
unstable cylinders of orbits which initially move towards the first bounce. In 
order to detect the first intersections of the two cylinders we adopt E : (x — 
Xb,Px = 0) as the surface of section, where Xb — 1.3118 is the coordinate of 
the first bounce of the orbits at p x = 0. In Fig. [7] we plot the points (zb,p z J 
of the sections of both cylinders in the first cross of E. The points A, B , C, 
D (contained in the sector (z,p z ) of E) characterize the transversal crossing 
of the cylinders. An detailed examination of the numerical points of the map 
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Figure 7: (first experiment) The first crossing of the stable cylinder (gray) 
and unstable cylinder (black) in the surface of section £ (at the first bounce) 
shown in the plane ( z,p z ), for E 0 = 7.9096. The four points A, B, C , 
D are the unique points of the first transversal crossing of the cylinders. 
These points define homoclinic orbits which are contained in both unstable 
and stable cylinders and are bi-asymptotic to the center manifold S 3 (E 0 ), 
constituting an invariant signature of chaos in the dynamics. Here E crit = 
7.909653935312942 
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shows indeed that all orbits arrive at the first bounce (xb ~ 1.3118, p x = 0) 
with coordinates Hi, — 1, p yb = 0, being a further verification of the accuracy 
of our numerical treatment. The points A, B , C, D in Fig. [^therefore define 
homoclinic orbits, namely, orbits which are in the intersection of the unstable 
and stable cylinders, and give an invariant characterization of chaos in the 
model. These homoclinic orbits are contained both in the stable cylinder and 
the unstable cylinder and are bi-asymptotic to the center manifold 5 b 0 C 
S 3 (E 0 ). 

The coordinates (y,p y , z,p z ) of the homoclinic points A, B , C, D are 
given approximately by 


A ~ (1,0,0.9997818,0.0330003), 

B ~ (1,0,1.0604379,0.0012771), 

C ~ (1,0,0.9991485,-0.0331563), 
D ~ (1,0,0.9445269,-0.0011262). 


Analogously in the second experiment we maintain the same values for 
the parameters Eq, k 2 and A together with the initial conditions (xo = x cr u , 
Px 0 — 0). However now fix the remaining initial conditions on the 4-dim 
invariant submanifod (ii) instead, namely, (y 0 = z 0 ,p yo = 3 p zo ). In fact 
we are restricting ourselves to a particular domain of initial conditions of 
the center manifold S' 3 which has the topology of S 1 and is defined by the 
Hamiltonian constraint 

H(x = x cr ,p x = 0,y = z,p y = 3 p z ) = 0. 

With the exact dynamics we generate one 2-dim stable and one 2-dim unsta¬ 
ble cylinders which initially move towards the first bounce. These cylinders 
are generated from a set of about 1, 500 orbits, with initial conditions taken in 
the above domain which actually correspond to a flow in the 4-dim invariant 
submanifold (ii) of the full 6-dim phase space. 

In Fig. |S]we plot the points ( Zb,p Zb ) of the sections of both cylinders in 
the first crossing of £. The four points E, F, G, H of the figure characterize 
the transversal crossing of the cylinders, defining homoclinic orbits which are 
in the intersection of the unstable and stable cylinders and are bi-asymptotic 
to the center manifold C S 3 (E 0 ). A complementary map in the sector 


31 
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Figure 8: (second experiment) Projection on the plane (z,p z ) of the first 
crossing of the stable cylinder (gray) and unstable cylinder (black) in the 
surface of section £ (at the first bounce). In this projection the four points E, 
F, G, H are the unique points of the first transversal crossing of the cylinders 
(cf. text). These points define homoclinic orbits which are contained in both 
unstable and stable cylinders and are bi-asymptotic to the center manifold 
S 3 (Eq), constituting an invariant signature of chaos in the dynamics of the 
model. 
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Figure 9: (second experiment) Projections on the complementary plane 
(■ y,p y ) showing the first crossing of the stable cylinder (gray) and unsta¬ 
ble cylinder (black) in the surface of section £ (at the first bounce). The 
four points E, F, G, H of figure El characterizing the transversal crossing of 
the cylinders are shown, constituting an invariant signature of chaos in the 
dynamics of the model. 
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(■ y,p y ) is given in Fig. El showing the first crossing of the cylinders by surface 
of section £ at the bounce. This map confirms the transversal crossings at 
the points E , F, G, H. The coordinates ( y,p y , z,p z ) of the homoclinic points 
E, F. G, H are given approximately by 

E ~ (1.0002743,0.0496506,1.0002743,0.0165502), 

F ~ (1.0290076,-0.0013846,1.0290076,-0.0004615), 

(58) 

G ~ (1.0002743,-0.0496506,1.0002743,-0.0165502), 

H ~ (0.9721427, -0.0009974, 0.9721427, -0.0003324). 

An detailed examination of the numerical points of the map shows that 
all orbits arrive at the first bounce (x — 1.3118, = 0) with coordinates 
(Vb = Zb) and ( p yb = 3p Zb ). This is also illustrated in Figs. flOl where the 
hrst crossing of the unstable and stable cylinders with the surface of section 
£ in the bounce are displayed. The section of both cylinders, projected on 
the sectors (y,z) and (p y ,p z ), he on the straight lines y = z and p y = 3 p z 
respectively, as expected. This is also a further verification of the accuracy 
of our numerical results. 

The sets (E7D and (1551) are two distinct numerical evidences of chaos in 
the dynamics, and constitute an invariant signature of chaos in the model. 
We must mention that the dynamics near homoclinic orbits is very complex 
associated with the presence of horseshoe structures[28], :27, j29j (22( (46]. The 
coordinates of the homoclinic points (157|) and (I58|) satisfy the constraint 
(I56[) . implying that they are contained in the transversal intersection of two 
S 3 at the bounce, namely, a S 2 . This fact indicates that the chaotic saddle- 
connected with the structure of homoclinic orbits bi-asymptotic to the center 
manifold S 3 - is a Cantor-type set having the compact support S 2 [35]. 

8 Conclusions and final comments 

In this paper we examined the dynamics of a Bianchi IX model, with three 
scale factors, sourced by a pressureless perfect fluid in the framework of 
bouncing Braneworld cosmology. Assuming a timelike extra dimension and 
a 5-D de Sitter bulk, the modified Einstein’s held equations on the 4-dim 
Lorentzian brane furnish a dynamics with correction terms that avoid the 
singularity and implement nonsingular bounces in the early phase of the 
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Py 

Figure 10: The first crossing of the stable and unstable cylinders with the 
surface of section £, projected on the phase space sections ( y , z) and (p y ,p z )- 
As expected the sections of both cylinders coincide on the straight lines y = z 
and p y = 3 p z . The points corresponding to the transversal crossing of the 
cylinders E, F, G and H are plotted for reference. 
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universe. In terms of metric functions the modified Einstein’s equations 
have a first integral that can be expressed as a Hamiltonian constraint in 
a 6-dim phase space, yielding a three degrees of freedom dynamical system 
which governs the motion in phase space. Due to an effective cosmological 
constant on the brane the phase space presents two critical points in a finite 
region of the phase space, a center-center-center and a saddle-center-center, 
plus two critical points at infinity corresponding to the de Sitter solution. 
Together with a 2-dim invariant plane of the dynamics the critical points 
allow to organize the dynamics of the phase space. 

We examine the structure of the dynamics in a linearized neighborhood 
of the saddle-center-center. We identify constants of motion associated with 
the saddle sector, which allow to define the linear stable V$ and unstable 
Vjj manifolds. We also identify constants of motion connected to the center- 
center sector, which define the center manifold of linearized unstable periodic 
orbits and has the topology of S 3 . In the linear domain the direct product 
Vs x S 3 and Vjj x S 3 define the structure of stable and unstable cylinders 
which constitute boundaries in the 5-dim energy surface of the dynamics. 

The nonlinear extension of the center manifold of unstable periodic orbits 
is parametrized by the constant of motion E 0 (E 0 < E cr ) with the topology 
of S 3 maintained. As one decreases the parameter Eq the nonlinearity of the 
center manifold increases, with a corresponding increasing of the dynamical 
instability as shown in our numerical simulations. The extension of the 4-dim 
stable Ws = Vs x S 3 and unstable Wjj = Vjy x S' 3 cylinders away from the 
neighborhood of the center manifold have the structure of 4-dim spherical 
cylinders with the topology R x S 3 . 

By a proper canonical transformation we are able to separate the three 
degrees of freedom of the dynamics into one degree - connected with the 
expansion and/or contraction of the scales of the model - isolated from the 
other two related to pure rotational degrees of freedom associated with the 
center manifold S 3 . By expanding the Hamiltonian constraint and Hamil¬ 
ton’s equations in these coordinates we show that the typical dynamical flow 
is an oscillatory mode about the orbits of the invariant plane. In particular 
the stable Ws and unstable Wu cylinders are composed of oscillatory orbits 
about the separatrices which emerge from the saddle-center-center critical 
point and guide the cylinders. These cylinders have the same energy Eq of the 
center manifold and coalesce to it as t —> ±oo. As these spherical cylinders 
are 4-dim surfaces they separate the 5-dim energy surface in two dynamically 
disconnected pieces. This fact is a fundamental feature of the dynamics for 


36 


characterization of chaos in the case of an eventual transversal intersection 
oiW s and Ws- As the separatrix which divides regions I and If in the invari¬ 
ant plane makes an homoclinic connection to the saddle-center-center critical 
point, this fact necessarily leads to the transversal crossings of Ws and Wu- 
The transversal crossing of the cylinders consists of homoclinic orbits which 
are contained both in the stable and the unstable cylinder and are biasymp- 
totic to the center manifold S 3 . The presence of a homoclinic orbit in the 
dynamics is an invariant signature of chaos in the model[28l [27| [23, [31]. The 
homoclinic intersection manifold has the topology of R x S 2 and consists of 
all homoclinic orbits biasymptotic to the center manifold defining a chaotic 
saddle^Mh associated with S 3 . 

The first transversal crossings of the stable Ws and unstable Wjj cylin¬ 
ders are shown numerically in two distinct experiments. For the sake of 
computational simplicity we restricted ourselves to cylinders generated from 
initial conditions taken on the center manifolds of the two 4-dim invariant 
manifolds of the dynamics defined respectively by (y = 1 ,p y = 0), and 
(y = z,p y = 3 p z ). We adopted the surface of section £ at the bounce defined 
by (xb,p x = 0) where Xb is the scale factor of the bounce for the orbits. By 
performing the evolution of orbits via the 6-dim exact dynamics we generate 
one 2-dim stable and one 2-dim unstable cylinders of orbits, and detected 
their transversal intersection corresponding to four homoclinic points in the 
first crossing of £ by the cylinders, for both experiments. These points define 
orbits which are homoclinic to the center manifold ^Se 0 C S 3 (E 0 ). 

In all our numerical simulations we used the 6-dim exact dynamics, in 
accordance with ([23]) . and the error in the Hamiltonian constraint (T5Tj) is 
checked to remain < 10~ 13 for all t. 

We finally compare some features of the dynamics, namely the oscillatory 
approach to the bounce and the chaotic behavior of the dynamics, with anal¬ 
ogous features present in the BKL conjecture in general relativity. First we 
note that in both models the oscillatory approach to the bounce/singularity 
is a key feature of the dynamics. In the general Bianchi IX model discussed 
here the three degrees of freedom of the dynamics are separated into one 
degree (connected with the expansion and/or contraction of the scales of the 
model) plus pure rotational degrees of freedom associated with the center 
manifold S 3 . The typical dynamical flow is an oscillatory mode about the 
orbits of the 2-dim invariant plane; in particular from the center manifold 
it emerges the stable and unstable 4-dim cylinders of oscillatory orbits that 
are guided towards the bounce by the separatrix in the invariant plane. In 
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the limit of k 2 —y 0 (close to the general relativity dynamics) the motion on, 
or about the unstable cylinder approximates the oscillatory BKL motion up 
to a scale x 3 > n 2 . As one can make n 2 as small as wanted, a long oscil¬ 
latory approach towards a neighborhood of x = 0 can be developed, with a 
behavior analogous to one of the Kasner eras of the BKL model. However 
k 2 cannot be made equal to zero as this would correspond to a change of 
topology of the phase space. The same considerations would apply for the 
case of a mixmaster universe, with a nonvanishing cosmological constant and 
Eq = 0, in general relativity. Second, the chaos in the present model has a 
homoclinic origin, resulting from the homoclinic transversal intersections of 
the stable and unstable 4-dim cylinders emerging from the center manifold 
S' 3 . In contrast the chaotic behavior in the BKL dynamics appears in a map 
that connects the length of the succeeding Kasner eras in the approach to the 
singularity of general relativity, for which we have no counterpart. Never¬ 
theless, considering the general relativity limit, we have topological evidence 
that the 4-dim cylinders - emerging from the center manifold S 3 and guided 
by the separatrix connecting the saddle-center-center to the singularity - 
should intersect and generate a homoclinic orbit from this intersection. 

In a future work we intend to examine the transversal intersection of the 
spherical cylinders R x S 3 in the full 6-dim phase space. We also intend to 
examine the chaotic exit to the final accelerated de Sitter stage for initial 
condition sets (corresponding to initially expanding universes) taken in a 
small neighborhood about the separatrix S approaching the saddle-center- 
center for t > 0. As in [13], we expect these sets to have fractal basin 
boundaries connected to the code recollapse/escape leading to a chaotic exit 
to the de Sitter accelerated phase. We also expect to observe the draining of 
initial condition basins from recollapse to escape behavior, as time increases. 
For t —> oo only the homoclinic intersection manifold is expected to remain 
in recurrent oscillatory motion. 
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